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Abstract: We examine squark-gluino loop effects on the process 6 — > 57 in minimal 
supersymmetry with general flavor mixing in the squark sector. In the regime of heavy 
squarks and gluino, we derive analytic expressions for the beyond-LO corrections to the 
Wilson coefficients and find them to be often large, especially at large tan/3 and ^ > 0. 
The ensuing ranges of values of the Wilson coefficients are typically smaller than in the 
LO approximation, and sometimes even change sign. This has the effect of often reducing, 
relative to the LO, the magnitude of supersymmetric contributions to BR{B — > Xg'^). 
This "focusing effect" is caused by contributions from: (i) an RG evolution of the Wilson 
coefficients; (ii) a correction to the LO chargino contribution to the Wilson coefficients, 
which can considerably reduce the LO gluino contribution. This partial cancellation of 
the two contributions takes place only in the case of general flavor mixing. As a result, 
stringent lower bounds on the mass scale of superpartners, which apply in the case of 
minimal flavor violation, can be substantially reduced for even small departures from the 
scenario. The often disfavored case of ^ < can also become allowed for MgusY as small 
as ~ 200 GeV, compared to > 500 GeV at LO and over 2 TeV in the case of minimal 
flavor violation. Limits on the allowed amount of flavor mixing among the 2nd and 3rd 
generation down-type squarks are also typically considerably weakened. The input CKM 
matrix element can be larger than the experimental value by a factor of ten, or can 
be as small as zero. 
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1. Introduction 

It has long been recognized that the inclusive process B Xg'j plays a prominent role 
in testing "new physics" beyond the Standard Model (SM) Q. In the SM the leading 
contribution comes from a virtual W^-top (charm, up) quark loop, which is accompanied 
by a real photon with GIM suppression. New physics effects, which are also loop induced, 
can be of the same order if a mass scale associated with the new physics is not much larger 
than mw- This in particular is the case with softly broken low-energy supersymmetry 
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(SUSY) where the scale of SUSY breaking MsusY is expected to be < 0(1 TeV) on the 
grounds of naturalness. 

Over the last decade, experimental precision has improved considerably, and the world 
average of the measured branching ratio has now been determined with an uncertainty of 
some 10%, 

BR{B ^ Xs-i)^^pt = (3-34 ± 0.38) x 10"^ (1.1) 

where results from four experiments [^, including a recent one from BaBar, have been 
taken into account Q. 

Theoretical calculations in the SM have been performed in several steps [Q, ^, including 
the full next-to-leading order (NLO) QCD correction which has recently been completed 
in and have reached as similar precision^ 

BR{B Xs-f)sM = (3-70 ± 0.30) x 10~^ (1.2) 

Given the approximate agreement of the SM prediction ( |1.2| ) with experiment ( |1.1| ), any 
possible "new physics" effects must now be confined to the remaining, relatively narrow, 
window of uncertainty. This is normally expected to impose severe constraint on extensions 
of the SM. For example, in the two Higgs doublet model (2HDM) extension of the SM, where 
full NLO QCD corrections have also been completed [0, ^, this allows one to derive rather 
stringent lower limit on the mass of the charged Higgs. It is worth noting, however, that 
the limit becomes sizeably weaker when NLO QCD corrections are included, as compared 
with the leading order (LO) approximation This demonstrates the importance of NLO 
QCD corrections [|^, ^] and NLO mass renormalization effects 

Supersymmetric contributions to BR[B ^s7) have not yet been calculated with a 



similar level of accuracy, although, since the first detailed LO analysis [10|, much important 
work has been done towards reaching this goal Q. On the other hand, one might argue 
that, given little room left for "new physics" contributions, it may actually not be essential 
to include NLO QCD and SUSY QCD corrections since all the relevant superpartners may 
anyway have to be heavy enough in order for the SUSY contributions to be suppressed. 
However, in this regime, the magnitude of SUSY contributions is often comparable with 
LO and NLO QCD corrections. Furthermore, in some cases separate SUSY effects may 
actually be quite large, but may approximately cancel each other, thus allowing for lower 
-^SUSY- 

More importantly, new effects are known to exist beyond LO, and these may signifi- 
cantly modify SUSY contributions which in the LO approximation appear to be large. One 
such important new effect is induced by squark-gluino loop corrections to quark masses 
and couplings. For example, the mass correction to the bottom quark can be of order 



40% in the regime of large tan /? |11|. The same mechanism has also a strong impact on 



effective quark couplings to gauge and Higgs bosons and to superpartners O, 13, 14, 15 



Its effect on 6 — > S7 has been explored [|14| , 15 1 in the framework of minimal flavor violation 



(MFV) [|l6[ (where one assumes that the CKM matrix is a unique source of flavor mixing 



^Somewhat difFerent ranges are quoted in the literature but they agree within their 1 a errors. 
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among both quarks and squarks), and shown to be sizeable at large tan/3 and for large 
ratio of MsusY/?7^w• 

The MFV scenario has some appealing theoretical features and is also partially moti- 
vated by stringent constraints from the — system [17, 2^, 18] on the allowed mixings 



among the squarks of the first two generations. On the other hand, it is hard to believe 
that MFV is strictly obeyed at the electroweak scale. It is usually assumed to hold at some 
high scale A, like the grand unification scale, in an attempt to relate soft SUSY breaking 
terms to the structure of the Yukawa sector. However, the assumptions of MFV are not 
RG-invariant^ and sizeable flavor-violating terms of order In (A^/m^) / (47r)'^ are gener- 
ated by the running from the scale A to the electroweak scale. Further violations are also 
induced by threshold effects at the SUSY breaking scale MsusY ~ 1 TeV. It is therefore 
sensible to assume that MFV can be regarded as, at best, only an approximate scenario at 
the electroweak scale. 

A broader framework is that of general flavor mixing (GFM) among squarks and/or 
among sleptons. Flavor mixings of this type are not protected by any known symmetry, 
and can in principle be large. As a result, without any additional ansatz, otherwise allowed 
contributions to, e.g., the processes such as ^ — > 67 or mixing can exceed exper- 



imental bounds by as much as several orders of magnitude |21|. In this case, experimental 
limits become an efficient tool in severely restricting possible mixings among scalar super- 
partners, especially between the 1st and 2nd generation, while still leaving considerable 
room for departures from the MFV scheme. 

In the case of 6 — > S7, going beyond the MFV framework has particularly important 
implications. While in the MFV case only the chargino and up-squark loop exchange 
contributes to the process, in going beyond MFV flavor-violating gluino (and, to a lesser 
extent, neutralino) and down-squark loop exchange diagrams start playing a substantial 
role. The importance of these additional diagrams has long been recognized p^], in part 
because of the gluino exchange enhancement by the strong coupling constant. Following 



the earlier work of |23| , 17|, it has more recently been carefully studied in the LO approx- 
imation [^], and also in some special cases beyond LO, e.g. in |2^]. (For a clear and 
systematic discussion of SUSY contributions in the case of GFM see ref. ||2^ .) 

Given the fact that beyond-LO effects in 6 — > 57 have been shown to be sizeable 
in MFV and that by going beyond MFV new important loop contributions appear, it is 
worthwhile to consider the role of beyond-LO effects in the framework of GFM among 
squarks, and to examine their impact on bounds on super symmetric masses and other 
parameters. In particular, one may ask whether small departures from MFV assumptions 
lead to only small perturbations on the various bounds on SUSY derived in the case of 
MFV. We will show that this is not the case. 

In the previous paper [^] we have presented our first results showing a strong impact 
of the squark-gluino loop corrections on 6 — > 57 in the framework of GFM. We have 
pointed out that such beyond-LO effects have an important effect of generally reducing 
the magnitude of SUSY contribution to BR{B Xg'^), especially at large tan/3 and for 



^For a comprehensive discussion and model-independent approach see ]l£ 
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// > 0. This effect of "focusing" on the SM value is already present at some level in the 
case of MFV but becomes strongly enhanced when GFM is allowed. This is because in this 
case large beyond-LO contributions from the chargino and the gluino can partially cancel 
each other. Furthermore, gluino-squark corrections typically reduce SUSY contributions to 
Wilson coefficients relative to the LO approximation. This effect comes in addition to 
the known effect of reducing the magnitude of the Wilson coefficients due to a resummation 
of QCD logarithms in RG-running between MsusY and mw [14|. In addition, in the 



framework of GFM, large contributions to the chirality-conjugate Wilson coefficients Cj g 
due to chargino exchange are induced which can partially cancel the contribution from 
gluino exchange. 

In this study, we present analytic expressions for dominant beyond-LO corrections to 
the Wilson coefficients for the process 6 ^ 57 in the Minimal Supersymmetric Standard 
Model (MSSM) with GFM in the squark sector. In addition to including the NLO con- 
tributions from the SM and the 2IIDM and the full NLO QCD corrections from gluon 
exchange, we compute dominant NLO-level SUSY QCD effects to relevant SUSY contri- 
butions, which are enhanced by large-tan (5 factors. While the full two-loop calculation 
still remains to be performed (with the first step already undertaken in [^), our goal is to 
calculate the terms which are likely to be dominant beyond LO due to being enhanced at 
large tan (3 and at large Msusy/"Iw- We will work in the regime where all the colored su- 
perpartners are heavier than all the other states. We generalize the analyses of [28, 14, 15| 



to the case of GFM. In addition, we include in our analysis an NLO anomalous dimension, 
an NLO QCD matching condition and a resummation of tan /^-enhanced radiative correc- 
tions which all further improve the accuracy of our analysis. We use the squark 
mass eigenstate formalism which is more appropriate when squark off-diagonal terms are 
not assumed to be small. We thus do not use the mass insertion approximation, similarly 
as in ref. p^] . 

Our main results can be summarized as follows. In the framework of GFM, dominant 
tan /3-enhanced beyond-LO effects due to squark-gluino loop corrections are quite large 
and have a strong impact on the Wilson coefficients. Their main effect is that of: (i) 
reducing, relative to LO, the size of the Wilson coefficients C/^S) and (ii) of generating 
additional contributions to the Wilson coefficients C7 g which can be of similar size. As a 
result, the total values of Cy'g can significantly change. For fi > these are usually have 
a smaller value, but not necessarily magnitude, relative to the LO. In fact, they can even 
change sign. (For fi < the pattern is less clear but the effect often remains sizeable.) As a 
result, the supersymmetric contribution to BR{B Xs'f) can often change dramatically, 
even for small mixings among down-type squarks. 

Our analysis demonstrates that various bounds on SUSY parameters obtained in the 
context of the MFV framework at the same level of accuracy can be highly unstable against 
even even small departures from it. An additional effect is that of strongly relaxing upper 
limits on the allowed amount of mixing in the 2nd-3rd generation squark sector relative 
to the LO case. Another important implication of our analysis is that the elements of the 
input CKM matrix can be larger than the experimentally measured values by a factor as 
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large as ten (which is significantly more than some ~ 40% obtained in [35| with a specific 
boundary condition at the GUT scale), or can be even zero, in which case the measured 
values of the CKM matrix would have a purely radiative origin [^] . 

The impact of large gluino-squark loops on Wilson coefficients, considered here for the 
process of 6 — > 57, appears to be quite generic and is likely to play an important role in 
other rare processes involving the bottom quark pO| ]. 

The paper is organized as follows. In sec. |2| and in appendices |A|-p| we provide a 
detailed discussion of the effect of the squark-gluino loop corrections on quark masses and 
couplings and on the squark sector. In sec. ^ we outline our procedure of computing the 
Wilson coefficients, and in sec. ^ we present their analytic expressions. Numerical effects 
are then presented in sec. |5| and our conclusions are summarized in sec. ^. Several useful 
formulae are collected in appendices 0^ for completeness. 

2. Effective Quark Mass in SUSY with General Flavor Mixing 

We will first define an effective, supersymmetric, softly-broken model with general flavor 
mixing in the squark sector as the framework for our analysis. In particular, we will concen- 
trate on the effect of the squark-gluino loop corrections to down-type quark mass matrix 
on other properties of the model, in particular on the squark sector, effective couplings and 
the CKM matrix. Unless otherwise stated, we will work in the MS'-scheme. 

The quark mass terms in an effective Lagrangian in the super-CKM basis read 

-C'^'''' = Sf^m^a^dL + d)j^6mddL + u^j^m^^^UL + u^j^SmM + h.c, (2.1) 

where dL = (dLi), dn = {dju), are the down-quark fields, and = [uli) and ur = {um) 
are up-type quark fields, with m-^^^ = {m'^^^ = (y^^d^ij) (^'J ~ 1)2,3) denoting their 
respective uncorrected, or "bare", 3x3 mass matrices. These are related to the respective 
input, or "bare", Yukawa coupling matrices in the usual way m^^^^ = Vd,uY^^^ where 

The mass corrections Smd,u arise from squark-gluino one-loop contributions.^ Of 



particular interest is the mass correction Smd since it is enhanced at large tan/3 [11| 
and/or in the case when left-right soft terms are large, 

Srud = Vd tan p Sy}^^ + Sm^f^K (2.2) 

Note that the latter are in general independent of the Yukawa coupling and can also be 
sizeable.^ An explicit expression for Sma and a self-consistent procedure for computing 
it will be given below. Analogous corrections in the up-quark sector are not enhanced by 
large tan /3-factors, although in general there could be a substantial effect from left -right 
soft terms. 



■^For simplicity we neglect all weak and Yukawa contributions as sub-dominant |^]. They can be found 
in ||. 

^In this paper, we adopt a purely phenomenological approach at the electroweak scale and do not assume 
any thoeretical constraints comming from, e.g., renormalization group evolution above the electroweak scale 
and/or vacuum stability conditions [p2[. 
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The effective (pfiysical) quark mass matrices mci,u in the super~CKM basis are given 

by 



ma = m^a^ + (^^^ Sma, (2.3) 

= mW + (^) 6m^. (2.4) 

Note tfiat we iiave expiicitely extracted the factor a/An in the expressions above in order 
to stress the order of the leading SUSY QCD (SQCD) one-loop correction. Subdominant 
effects, e.g., due to chargino-stop loops considered in can be added linearly. 

The effective (physical) quark mass matrices Tnd,u in the super-CKM basis are (by 
definition) diagonal: 

rud = diag {md^5ij) = diag {nid, ms,mb) , (2.5) 



ruu = diag {rriu^dij) = diag {rriu, rUc, mt) , (2.6) 

where md,s,b and mu,c,t denote the physical masses of the quarks. In contrast, in this basis 
the matrices will in general be non-diagonal, and so will 6md,u- 

The CKM matrix is defined as usual^ 

K = V^,VI. (2.7) 

The elements of K are determined from experiment, except for a small loop correction which 
can be found in appendix |^. The unitary matrices Vd,u^ ^ rotate the interaction basis' 
quark fields d°^^, u°j^-, u°j^- to the mass eigenstates dn, dRi, uli, um {i = 1,2,3) in the 
physical super-CKM basis. For example, dn = Vd^^jdl-, um = VuRijU°^j = 1,2,3). 

Also Tn^^ = "^d'^" ' '^here Tn^^" stands for the down-quark "bare" mass matrix 
in the interaction basis, and analogously for the up-quarks. 

For our purpose, it will be convenient to further introduce a "bare" super-CKM basis. ^ 
It is the basis in which the "bare" mass matrices "fn^lj^ are diagonal. In the absence of 

the mass corrections Sfn^^, the "bare" and the physical mass matrices would of course 
coincide and their eigenvalues would be equal to the physical quark masses which we will 
treat as input. In other words, there would be no need to distinguish between the "bare" 
and the physical super-CKM basis. (In fact, this level of accuracy is sufficient to compute 
supersymmetric contributions to rare decays in the LO approximation.) Note however, 
that once the mass corrections are switched on, it will be the "bare" Yukawa couplings, 

^We will generally follow the conventions and notation of Gunion and Haber unless otherwise stated. 

®To be precise, one should use different symbols to denote the "bare", or input, quantities, which may 
be expressed in different bases, in particular in the interaction basis and the super-CKM basis, and the 
various quantities, "bare" or loop-corrected, appearing in the "bare" super-CKM basis. For the sake of 
simplicity, we will use just one superscript "(0)" but will define all the relevant terms as they appear. 
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which initiahy appear as free parameters in the theory, that wih have to be adjusted to 
compensate for the loop corrections. 

The "bare" CKM matrix is defined as = VS^^v}"^ \ where v}^?, denote the 

unitary matrices that would transform the "bare" down- and up-quark mass matrices in 
the interaction basis into a diagonal form, in analogy with the case of the corrected super- 
CKM basis. As we will see later, in the presence of the one-loop mass corrections, the 
elements of and Kcb (the latter being determined by experiment) can show large 

discrepancies, which can reach a factor of ten, or can be reduced down to zero. 

Turning next to the squark sector, in the super-CKM basis spanned by the fields 
d-L = (dLi), dii = (dijj {i = 1,2,3), and ul = {uLi), ur = {ur^), the mass terms in the 



Lagrangian read 




into 3x3 block sub-matrices as follows 



The 6x6 down-squark mass matrix ~ and the up-squark matrix A4~ are decomposed 



^2 _ I "^XlL + Fd,LL + Dd,LL rnl LR + Fd,LR 

' ^^d,LR + ^d,LR) rnl 

,RR + ^d,RR + Dfi,RR 



'^u,LL + Fu,LL + Du,LL 



^l,LR + Fu,LR 



^u,LR + ^u,LR 



^u,RR + Fu,RR + Du,Rii 



(2.10) 



The terms appearing in ( |2.9D -( 2.10 ) are related to their more familiar counterparts in the 
interaction basis by the same unitary transformations Vd,u^ that appear in the quark 
sector: squark fields transform in the same way as the respective quark fields from the 



interaction basis <i£ 



Hi 



u^J to 



(i = 1,2,3), and m£ = (n°J, Wj^ 

the super-CKM basis which is spanned by <ii,,H ai^d ul^h. For example, dL = Vd^^d^ and 
Ul = Vuj^u^, etc. If one writes the Lagrangian for the soft SUSY breaking terms in the 
interaction basis as 



^o, mass 
q, soft 



dVmldl + drJmld<}^ + 



dV {^dAl' 



d^ + h.c. 



+ 



Ur ' [VuA 



(2.11) 



then the soft mass terms appearing in the mass matrix A4'i are given as fohows 

d 

mlLL = y<i,mlvl, (2.12) 

mlaR = Vd^m^oyL (2-13) 

mlLR = y<iA^dA*^)yL (2-14) 
and similarly for the up sector. (In particular, j^j^ = ^h'^^V^^). The matrices 

iT^du LL RR hermitian and are in general non-diagonal, as are their 
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counterpart soft mass parameters "ttIq, m/^ and in ( 2.11| ). Note that rn^ j^j^ and 
are related by SU{2)l invariance. Since in the interaction basis they are both 
equal to ?tIq, in the super-CKM basis one finds = K^m^ ll^- 

The trilinear 3x3 matrices Ad^u that appear in the LR soft terms are in general 
non-hermitian. Note that we do not assume -^d^ti to be necessarily proportional to their 
respective Yukawa couplings. 

The F-terms are given by 

Fa,LL = m'^d^^m''d\ (2-15) 
Fd,RR = m'^a^m'^a^\ (2.16) 
Fd,LR= -fJ^ tan /3m'^^^\ (2.17) 

and similarly for the up squarks (except Fu,lr = — fJ. cot P mj^^^ ) . Note it is the "bare" 
mass matrices that appear in the F-terms. This is because these are derived from 

the superpotential which contains "bare" Yukawa couplings and other "bare" parameters. 
For this reason, in the super-CKM basis the F-terms are not diagonal, which will lead to 
important effects. On the other hand, the D-tevms remain diagonal in flavor space. They 
read 

Dd,LL = mz cos 2/3 [Tgd - Qd sin^ dw] 1, (2.18) 
Dd,RR = {mz cos 2/3 Qa sin^ dw) 1, (2.19) 

where T^d = — 1/2 and Qd = — 1/3 (and analogously for the up-sector, with T^u = 1/2 and 
Qd = 2/3). 



matrices which are represented by two 6x3 sub-matrices TdL and TdR {^uL and T^r 



The mass matrices ~ (-^D in the super-CKM basis can be diagonalized by unitary 



irdL,rdR)Ml{ J^f^ J =diag(m|,...,m|J (2.20) 




and 

(r.L,r^fi)M? I >^ I =diag{ml,...,ml), (2.21) 



where m^^ and 111^1 {I = 1,...,6) are the physical masses of the down-, d = ydj 
and up-type, u = (uj), squark mass eigenstates, respectively. These are related to the 
corresponding fields dL^n, ul^h in the super-CKM basis by 

d = TdL,RdL,R, u = TulUl,r- (2.22) 

We will also need squark mass eigenvalues and diagonalizing matrices in the limit of 
vanishing mass corrections Smd,u- Squark mass matrices in this limit will be denoted by 

M^~^ ^ = Ml {Srud ^ 0) , (2.23) 
M^^^^ = Ml{Sm^^O). (2.24) 
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Since in this limit the "bare" and the physical super-CKM bases coincide, the F-terms 
are now diagonal and are given by F^^^j^ = F^^^j^ = diag (m^, m^, m^) and F^^j^ = 
—fi tan P diag {md,ms,mb) (and analogously for the up-squarks). This is because the 
"bare" mass matrices Tn^^, that appear in the F-terms, are in this limit diagonal, with 
the physical quark masses along the diagonal. 

We will further introduce the masses of the down- and up-type squarks, respectively, 
in the limit of vanishing mass corrections Smd,u 

mf^ = {5md ^ 0) , mi] = m^, {dm^ ^0) (I = 1, . . . , 6), (2.25) 

as well as the diagonalizing matrices r^°^ = ^dL {^i^d 0) and ^^^R — ^dR {^f^d 0)> 
and analogously for the up-squarks. 

We will now apply the above formalism to describe our procedure of computing the 
one-loop mass correction 6md ■ The relevant diagram is given by fig. [l| for which one 
obtains 

6 

{Sm,\^ = 2C2{3)J2inR)H {rdL)j,mgBo{ml,ml), (2.26) 
1=1 

where the quartic Casimir operator for SU{3) is C2(3) = |, rrig denotes the mass of the 
gluino and Bq is one of the Passarino-Veltman functions which are collected in appendix 
(An analogous expression for the mass correction {5mu)^ can be obtained from ( 2.26| ) by 
simply replacing d ^ u and d ^ u.) 

In our analysis we take the physical masses m^j 2 a input. On the other hand, in 
computing the diagram of fig. |l|, one uses the elements of the "bare" mass matrix 
and of the squark mass matrix Al ~, which in turn are functions of and other "bare" 

d ^ 

parameters. The final result must be equal to the physical quark masses. This means that 
we need to employ an iterative procedure to simultaneously determine 'm^\ Sm^\m^^) 
and A4~(m^°^) to a desired level of accuracy. 

In the first iteration one neglects Sma and thus identifies with ma- One then 

computes Snid by applying ( ^.26 ) where the elements of the squark sector are determined 
by •^j'^^; aiid not by In other words, in evaluating ( p.26| ) one initially makes the 

substitutions mj^ —>■ rri~~^ and TdRL —>■ ^^rl- 



dj 



dLj — » KqQQQQQSL^ "* ^Ri 



Figure 1: Gluino-squark correction to the down-type quark mass matrix. Direction of the arrows 
indicates chirality of the quark field. 
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In the second step one uses Srud from the first step to compute m^^^ = — 
{as/An) SrUfi- The new values of the elements of 'm^^ are then used in ( p. 91) and ( 2.20 ) to 
determine the new eigenvalues m^^ = 1, • • • i 6) and the new elements of T^l and r^H- 

These are then used to compute Sma via ( |2.26| ) and in the third step, and so on. 

The procedure converges quite rapidly, although not as well at very large tan/3. 

The above iterative procedure amounts to resumming tan/3 terms to all orders in 
perturbation theory. The importance of resummation has been stressed in [15, 14 1 for 
extending the validity of perturbative calculation to the case of large tan /3. Below we will 
provide a simple formula for tan /3-resummation effect on the bottom mass in the limit of 
MFV. 

Squark-gluino loop corrections also modify various couplings involving gauge and Higgs 
bosons, as well as those of the chargino, neutralino and gluino fields. The effective couplings 
are summarized in the appendices ^ and |^. They are computed in the same self-consistent 
way as the quark and squark masses and mixings. 

We will now illustrate the above procedure of computing quark mass corrections in 
the limit of vanishing inter-generational mixings among quarks. Concentrating on the 
correction to the mass of the b-quark one obtains 



dnib = {5md)^2. = (^2(3) sin 26*gmg Bq \^m~,m~J - Bq (^m~,m~J , (2.27) 

where now the two sbottom mass states are bi = and 62 = (Iq, and they result from the 
mixing of only the interaction states bi and bn, with 9^ denoting the mixing angle, 

<2 



r F \-( ""^^^b sin6^\ 0^.. _ ^ l"^rf/36 /n 9o^ 

idL,-l-dR)- -a /i: ' sm26/j-^ l^-^o) 

/ \ — sm Ur cos th I mA — m= 

where = ("T'Ilr) + i^d.LR)^:^- After some simple steps one obtains 

' 36 ^ ' 33 

(m% 

where the function H2 is given by 



2 


2 


mA 




bi 




2 


2 











5mb = 2C2{3)^^^H2 \^,^\, (2.29) 



In the regime of moderate to large tan /3, and after neglecting j^j^ , one recovers a 
familiar expression |14, 15 1 

( 

where rn~^^ is the "bare" bottom quark mass and e?, is given by 



/ y(0) \ 

rrib = V2mw-^ cos/3 (1 + ebtan/3) = m^^' (1 + efetan/?) , (2.31) 



Svr nig 



2 


2 


m= 


mA 


bi 


62 


2 


2 


m~ 


m~ 


9 


9 



(2.32) 
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The resummation of large radiative corrections at large tan/3 |15| , 14 1 to all orders in 
perturbation theory is in our case achieved by employing the iterative procedure described 
above. In the case of the b-quark mass considered here, after the first iteration one obtains 



^(0) ^-|^st-j _ _ ^g^j-^ ^ where e[^^ denotes the correction ( |2.32| ) obtained in the first 



step of iteration, where one makes the replacement m-r mi^^ . If one denotes by e^'^^ 



1,2 bi 



the correction ( 2.32 ) obtained in the fc-th step of iteration then, after n steps one obtains 



m[°V*^) = mb[l- e[") tan/3 + e["^e["-') tan^ /? + ... + (-1)" e^^^e^""') . . . e[') tan" P 

(2.33) 

We will now discuss the case of MFV as a limit of the GFM scenario. While in the 
literature the notion of MFV is not unique, and often depends on a model, the scenario 
is usually defined as the one where, at some scale, the LL and RR soft mass matrices 
are proportional to the unit matrix and the LR ones are proportional to the respective 
Yukawa matrices. A convenient basis here is the super-CKM basis [^] where the Yukawa 
matrices are diagonal. In this case MFV corresponds to all off-diagonal entries in the 
soft mass matrices being zero. Thus it is the off-diagonal elements of fn^uLL^ ^duRR 
and rn'ciu lr which in this basis are responsible for flavor change This is because, in 
the super-CKM basis, the couplings of the neutral gauginos are flavor-diagonal while the 
mixing in the charged gaugino-quark-squark couplings are all set by the SM CKM matrix 
K. 

As mentioned above, in the absence of quark mass loop corrections, there is no need to 
distinguish between the "bare" and the physical super-CKM basis. However, once these 
are taken into account, an ambiguity arises. In one approach, one can insist that the 
LR soft mass matrices remain proportional to the "bare" Yukawa couplings, which now 
become modified with respect to their tree-level values. This view can be motivated by an 
underlying assumption that, at some high scale, the two sets of quantities are ultimately 
related. Alternatively, in a more phenomenological approach like this one, one can assume 
that the LR soft mass matrices do not change and thus remain proportional to the physical 
Yukawa couplings (which are the same as the "bare" tree-level Yukawa couplings) which 
in the super-CKM basis give the physical masses of the quarks. Here we choose the latter 
option. Fortunately, the difference between the two choices is numerically small. 

We thus define the MFV scenario as the one in which, in the physical super-CKM 
basis, one has 



^d,LL 0^ ^%RR ^ ^d,LR ^ ^d/Vd, 



(2.34) 



and analogously for the up-sector. Note that in the last relation of eq. ( 2.34 ) one also 
assumes that oc Y(i,u which in the case of GFM does not have to be the case. 

In order to parametrize departures from the MFV scenario, it is convenient to introduce 
the following dimensionless (3x3 matrix) parameters = {Sfi)-^, ^lr^ where 
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hj = 1,2,3, 



as well as 



^1l).= ' ; . , Kr).= \ ^ (2-35) 



4r).= . . , f4L).. = ^— — ^^==, (2.36) 



^lRR)..m,RR) .. \ (^Irr) .\^d,LL, 



and analogously for the up-sector. The 5's can be evaluated at any scale equal to or above 
the typical scale of the soft terms. We will compute them at the scale MgusY which we 
will assume to be the characteristic mass scale for the soft mass terms. Note that the 
definitions ( 2.35| )-( 2.36| ) remain the same with and without loop quark mass corrections 



and are therefore particularly convenient for comparing the effects of LO and dominant 
beyond-LO corrections to BR{B Xs^) considered here.'^ 



3. Outline of the Procedure 

Before plunging into technical details, we will first provide a general outline of our procedure 
for treating beyond-LO corrections to BR{B ^s7) in the framework of GFM. One starts 
with an effective Lagrangian [^] evaluated at some scale fiQ, 

L = ^KlKt, [C^{^^io)P^ + C'MPl] . (3.1) 

The Wilson coefficients Ci{^Q) and C[{iiq) associated with the operators Pi and their 
chirality-conjugate partners P[ play the role of effective coupling constants. The C'^s are 
obtained from Cj's by a chirality exchange L <^ R. 

At ^0 = /^w, the Wilson coefficient of the four-fermion operator C2 has a tree level 
contribution. Dominant loop contributions determine the coefficients C^^s of the magnetic 
and chromomagnetic operators, respectively, which are also most sensitive to "new physics". 
These coefficients mix with each other when they are next evolved by renormalization 
group equations (RGEs) from //w down to /Ub. Then BR{B Xg'y) is evaluated at /ib as 
described in Q, with updated numerical values for the "magic numbers" taken from |^]. 

3.1 SM and 2HDM Contributions 

In the SM+2HDM, in addition to the W^-loop, there is a charged Higgs boson contri- 
bution to which always adds constructively to the one from the SM. Both contribute 

^Alternatively, we can define the dimensionless parameters in "bare" super-CKM basis in view of the 
correlation between the soft masses and the "bare" Yukawa coupling, as stated above. The iterative 
procedure we have introduced works with this definition equally well. 
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to C/ as well, but this is suppressed by ms/mt, and can be ignored. In contrast, this is not 
necessarily the case in SUSY extensions of the SM. 

In the SM+2HDM case the sum in eq. (|3.1| ) extends to eight, while in SUSY in general 



several additional operators arise [24|. In our analysis we will henceforth restrict ourselves 
to the "SM basis" of operators but will consider below to what extent it is justified to 
neglect the additional operators. 

Corrections from hard gluons to the vertices involving the M^-boson, H~ and G~ 
(where G~ represents the unphysical scalar appearing in gauge) are evaluated at the 
scale /Ltw = In computing them we follow the NLO QCD calculation of the matching 
condition at /iw 0i |§| • 

3.2 SUSY Contributions 

In SUSY, additional loop contributions arise from diagrams involving the chargino, neu- 
tralino and gluino exchange, along with appropriate squarks. Furthermore, as has been 



shown in |12, 15, 14], NLO SUSY QCD corrections to one-loop diagrams involving the 
gluino field can be as important as those involving hard gluons. 

The diagrams involving a loop exchange of the wino and higgsino components of the 
charginos (the superpartners of the W~ and the H~) with the up-type squarks add to the 
SM part constructively or destructively, depending on the sign of fi. The chargino/up- 
type squark exchange is the only SUSY contribution if one imposes MFV assumptions. 
In the case of GFM, the diagrams involving the exchange of the gluino (and neutralino) 
and down-type squarks with flavor violating squark mass mixings are also allowed and, as 
discussed in sec. |2|, in general play a substantial role. 

In the softly-broken low energy SUSY model, with either MFV or GFM, an additional 
complication arises that is related to the fact that several new mass scales appear which 
are associated with the masses of the contributing SUSY partners. In this analysis we 
will assume that all color-carrying superpartners (gluinos and squarks) are considerably 
heavier (by at least a factor of a few) than other states 

A^SUSY ~ 0{mg,mg,mj^,mjj > /iw ~ 0{m'w,mt,mH±) > /ib ~ 0{mi,). (3.2) 

Such a hierarchy is often realized in low-energy effective models with a gravity-mediated 
SUSY breaking mechanism because of different RG evolution of soft mass parameters from 
the unification scale down. It thus appears justified to associate the heaviest superpartners 
with the SUSY soft-breaking scale MsusY- In particular, we assume both superpartners 
of the top quark to be heavy. Below the scale /^SUSY = -^SUSY we will then deal with an 
effective model with the gluino and squark fields decoupled. Their effect will nevertheless 
be present in modifying the couplings of the states below nsvsY- 

As stated above, for the states at /iw {W, G~ and H~), the contributions to the 
Wilson coefficients are evaluated at /^w The beyond-LO SUSY QCD corrections to these 
vertices, which are of order 0{as), are absorbed into the effective vertices involving Yukawa 
couplings by integrating out the gluino field using the effective coupling method which has 



been introduced in |28]. In our case, the effective couplings are computed at ususY and 
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then evolved down to ;Uw using RGEs for the SM QCD /3-functions with six quark flavors. 
For simplicity, in the effective couplings we neglect the running of the Yukawa and gauge 
couplings other than as- We have extended to the case of GFM in the squark sector the 
calculation of the effective vertices and of the complete NLO matching condition at //w 
which in has been done in the framework of MFV. 

The contributions to the Wilson coefficients from the chargino, gluino and neutralino 
vertices, including NLO QCD and SUSY QCD corrections, are first evaluated at /isusY 
and then evolved down to /uw with the NLO anomalous dimension corresponding to six 
quark flavors [^], and with NLO QCD matching condition [41| at /^SUSY within the SM 



operator basis. The leading SUSY QCD corrections beyond LO to the vertices of the above 
superpartner have been computed in refs. 14, ^] in the case of MFV. Here we extend 



the calculation to the case of GFM among squarks. 

In light of the assumed mass hierarchy ( p.2[ ), the chargino, gluino and neutralino 
vertices all involve some massive (~ Mgusy) fields. In this case expansion by external 
momentum is not justified and the effective coupling method, which has been applied 



above to the states at //W) cannot be used anymore |28|. In order to treat them properly, 
one would have to calculate full 2-loop diagrams involving heavy fields, which is beyond 
the scope of this study. This has not yet even been done in the less complicated case 
of MFV. Instead, in |14, ^] finite threshold corrections to the bottom quark Yukawa 
coupling, which are enhanced at large tan/3, have been considered in the case of MFV. 
Here we extend this method to the GFM scenario and compute tan /3-enhanced (and also 
A-term enhanced) finite threshold corrections to the matrix of the Yukawa couplings of 
the down-type quarks. (Analogous corrections to the up-type quarks are not enhanced at 
large tan/?, although can possibly be by the soft terms.) These beyond-LO corrections are 
likely to be dominant due to being enhanced by tan/? factors relative to other corrections. 
As mentioned above, we introduce NLO resummation of QCD logarithms to C/^g and 
Cj g assuming the hierarchy (^]^) . In computing the chargino, gluino and neutralino loop 
contributions we follow in implementing NLO QCD matching conditions. 

The full set of operators in the sum in ( |3.1D that arises in SUSY has been systematically 
analyzed in ref. |2^]. Box diagrams involving the gluino field generate scalar and tensor 
type four-quark operators, which are absent in the SM, at the matching scale /^suSY- As a 
result of RG evolution from /LtsuSY to fih, in contrast to the vector type operators in the SM, 
at 1-loop level these new operators mix with the magnetic and chromo-magnetic operators 
and the contribution due to the mixing from these new operators could be comparable to 
the initial loop contribution to the (chromo-) magnetic operators. 

This mixing effect appears, however, to be subdominant at least at LO.*^ In the presence 
of the new operators, the (chromo-) magnetic operators are classified as dimension 5 or 6 
depending on whether the chirality flip of the operator originates from the mass of the 



gluino in the loop or from that of the external quark [24|.^ The new operators mix only with 



We thank T. Hurth for providing us with this argument. 
^In our numerical calculations, we do not distinguish these two contributions and use the same operators 
and corresponding anomalous dimensions as for the SM contribution. This is justified when the new 
operators are neglected [E4|. 
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the dimension 6 operators and their contribution to the b ^ s'j ampUtude is suppressed 
by mb/m-g relative to that of the dimension 5 ones. Because the new operators always 
appear with the dimension 5 operators, which are generated by the same set of flavor 
mixing among squarks, their effect is subdominant. This suppression has been numerically 
confirmed in |24]. 

In the absence of the NLO anomalous dimensions, matrix elements, bremsstrahlung 
corrections and matching conditions for the complete set of operators, we are unable to 
evaluate their contribution at NLO at the same level of accuracy as for the SM set. How- 
ever, the above suppression mechanism is based on the chirality structure and we can see 
no obvious reason for any sizeable enhancement from the new operators at NLO which 
would overcome an additional suppression with respect to the already subdominant LO 
contribution. Based on the above argument, in this work we restrict ourselves to the SM 
set of operators.^'' 



4. Wilson coefficients 

In this section, we present expressions for the supersymmetric contributions to Cj\fi\Y) 
and Cs{fJ-w) obtained using the procedure outlined above. 

In the MSSM, the Wilson coefficients for magnetic and chromo-magnetic operators Cy 
and Cg, in addition to the W~-t and H~~t loops, receive contributions from the xr2~^iv,6; 
Xi,...,4-f^i,...,6 and g-di^,,,fi loops, 

C7,8M = S'^Cj^sM + <5^C7,8(/iw) + S^'Cr^sM, (4.1) 
where S = x~,X°,9, 



S-'C7,8M = 5^"C7,8(/uw) + C7,8(mw) + S^^Cj^sM- (4.2) 

Analogous expressions apply to C'^ § except for the cases explicitly mentioned below. 

The coefficients S^Cj^s, where X = W,H, including order a, NLO QCD and SUSY 
QCD corrections are evaluated at /iw 

5^C7,8(/iw) = S^ci'lM + ^^S^ci%w), (4.3) 

where 

<5^CS and d'^ci^l denote the LO and NLO contributions, respectively. 



^"At large tan/3, some of the box diagrams are enhanced by tan^ /3, instead of tan/3, as in the penguin 
diagrams. However, they always come with a suppression factor of (mi,/mg)^ and we assume that these 
contributions remain to be subdominant even at large tan f3 if the gluino is heavy. We also neglect a neutral 
Higgs mediated contribution to the new operators, which could be important at large tan/3 [^. 

^^Our notation for the Wilson coefficients follows that of ref. but is extended to the case of GFM. 
Note that in our case we do not assume them to be necessarily proportional to CKM matrix elements. 
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The H and W contributions to the Wilson coefficients at are matched with an 



effective theory with five quark flavors. Their LO contributions are given by |42, 28 1, 

,(1) fm'^im) 



WJ 



7,8 



m 



w 



1 ( mKm) 

3tan2/3^7,8 \^ 



(2) / rn'^im) 



7,8 



(4.4) 
(4.5) 



where mt{iJ,-w) is the MS running mass of top quark at /iw and mn denotes the mass of 

/-I r)\ _ 

the charged Higgs. The mass functions -F^g are given in appendix^. 

The NLO matching condition of the Wilson coefficients d^'^C^^l at //w reads 

^^'^CrftVw) = 6Y'"ci^lM + 5f ^cg(Mw), (4.6) 

where S^'^ci^l represents an NLO QCD correction from two-loop diagrams involving one 



gluon line and analogously 6~ 



Cj o stand for corrections from two-loop diagrams with a 



gluino line, which are integrated out at ^SUSY- Explicit expressions for sY'^'ci'^s 

are given 



in refs. 41 1 and collected in appendix ^. 



As outlined in sec. ^, in computing S^'^Cj^l we use the effective coupling method [28| 
and replace in the relevant LO diagrams the tree-level couplings for the W~ {G~) and H~ 
vertices by corresponding effective couplings. These are given in appendix ^ and are first 
computed in the regime of eq. (^.21) at //susY and next evolved down to /iw using RGEs 
with SM QCD /3-functions. The effect of the running of the Yukawa and gauge couplings 
other than is neglected as subdominant. One finally obtains 



"9 ^7,8 



L L 



Gmdiil)* ^Gusd3{0)\ p(l) 
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f~liiu.3U,3yyi\ t-i 

'-'L ^7,8 
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Hu3d2{<S)* ^Hu3d3{l) 
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Hu3d2{l)* ^Hu3d3{0) 
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(2) 
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7,8 



m^(/iw) 
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where mtil^'w) and rnhifiy^) are MS running top and bottom masses at //w Explicit 
expressions for the effective couplings (j^"^^^]^^'^^ g^^(^ ^Gu^dj{o,i) other effective 

couplings are given in appendix Analogous expressions for S^'^Cj g can be obtained 
by simply interchanging L and R in the formulae above. 

In light of eq. (p.2|), the LO and NLO supersymmetric contributions to the Wilson 
coefficients S'^Cy^s are first computed at fJ-susY, 



<5^C7,8(//SUSy) = <5^4°gVsUSY) + 



sMo)f.. \ I «^(/^susy) ,sM'^)r.. ^ 

<--7,8(/^SUSYj 



(4.9) 
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Next, they are evolved down to /xw with the QCD renormalization group equation. For 



this purpose we use the NLO anomalous dimension obtained in |39], assuming six quark 
flavors, and employ the NLO expression given in for matrix evolution 

(4.10) 



<5^C7(//w) \ _ f^, asim) j\ f .(0) f-n _ Qs(^susy) ,\ / S^Cjinsusv) 



- ' 7^21 — rj3 



where rj = as(^susY)/os(Mw)- After extracting the LO and NLO parts of Wilson coeffi- 
cients at both scales, eq. ( |4.10D reduces to a more familiar form in the case of six quark 
flavors, 

^^Cf (Mw) = v'^^6'c!f\^,susY) + ^ (^^ - S'C^s Hf^susY), (4.12) 
S^C^M = r?^5^cf (^5C/5y), (4.13) 

/584 5 76256 37 14296 2 3520 leX ,c^(o). 

+ 773 -I 7721 Co (ususy) 

V 27 3381 ' 621 ' 147 ' y ^ yn-bubv i 

+r,i<5^C«(^sf/5y) + \{^^- V^) 5'c'i\iJisusY), (4.14) 

.5^(1). ^ _ ^1x5/^(1)^ A , 73 / I ^2N „(0), 



5^C^^^(mw) = ^-^5^Cl'\t^susY) + y - j <5^Q '(^st/5y). (4.15) 

While we work within the SM basis of operators, SUSY loop induced 0[as) contri- 
butions to Ci {i = 1,...,6) at /xsusY also in principle mix with Cy^g during the course 
of evolution. However, the corresponding Pi {i = 1, ... ,6) conserve chirality and are not 
enhanced by tan/3. We have numerically verified that their typical values are 0(10"'^)- 
0(10"^) if squarks/gluino masses are not very close to //w- Their contributions to Cy^s are 
further suppressed by mixing factors. Because the purpose of this paper is to present large 
tan /3-enhanced beyond-leading-order effects in BR{B Xg^), for simplicity we neglect 
these small contributions compared to the 0(10"^) SM ones. 

The NLO super symmetric contributions to the Wilson coefficients at /isusY reads 

5^C«(//susy) = S'gCi%susY) + 6lci%svsY) {S = x",X°,9), (4.16) 

where 6gC!^^l and S-C^^'^ denote respective gluon and gluino NLO corrections, analogously 

to the W and H contributions at //w- Explicit expressions for SgC^^g can be found in 
ref. |4^] and appendix ^. 

In computing gluino corrections to super symmetric contributions 5~Cy^g, the effective 
vertex method does not work, as explained in sec. ^. Instead of calculating the full two- 
loop diagrams, we take into account tan /3-enhanced corrections to the Yukawa couplings. 



^^In our numerical analysis, we actually use eq. (4.10) and and thus keep 0{a^) and higher order correc- 
tions in Cr^sifJ^w) that are induced by the evolution. 
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Figure 2: Leading diagrams that contribute to (left) and S^C'^^'l (right). The photon (or 

gluon) hne is attached in every possible manner. 



as well as LR mixings due to terms involving the trilinear soft terms A(i,u which in principle 
can also be sizeable. In the absence of full 2-loop calculation of such effective Yukawa 
vertices contain only the dominant corrections. (In fig. ^ and |3|, we show the Feynman 
diagrams that contribute the leading corrections to C73.) We indicate the approximation 
by replacing 5~Cj \ — s- S-C'^^ . We obtain 
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+ A[7.8r' 

(4.18) 



4vr 9 7„ 



where 



2,glKlKtb ^ 



rrig 



"T-6(/^SUSy) 



5 



[7,8] f^d.1 



(my)2 



(X = ui.di, Y = XayXr^g), 



(4.19) 



(4.20) 



and h\^'^q and AI^'^I are given in appendix Here ^^(//susy) denotes the running bot- 
tom mass at /isusY- The couplings {GdR,L)iav i^dR,L)iri and {GdR,L)ij correspond to 
the chargino, neutralino and gluino vertices including gluino corrections and are given in 
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Figure 3: Leading diagrams which contribute to S~ g (upper) and S-C!^^ (lower). 



appendix Note that in the mass functions we use physical squark masses m and niuj , 
where / = 1, . . . , 6. 

The LO contributions at //susY S^c!fl {S = x~jX^:9) can be obtained from the 
above expressions by replacing m^^ — > m^p and rriuj — > m-'^ above and by replacing the 
effective chargino, neutralino and gluino couplings with their tree-level values, as discussed 
in appendix Expressions for beyond-LO contributions 6^0^^^ can then be obtained by 

subtracting from the full expressions for (^5^0^^^ + j^^^Cj^^^ in eqs. ( [4.17| )-( p:.19| ). 

Finally, expressions for S^C'^g and 6^C'^^^ can be obtained by interchanging the indices 
L and R in the above formulae. 



5. Results 

We will now demonstrate the effect of beyond-LO corrections to the Wilson coefficients 
derived above and to BR{B ^s7) with some representative examples. We will work 
in the super-CKM basis and in deriving our numerical examples will use the following 
parametrization: {^^11)^ ~ "^f'^M; i^dLL)ij ~ ^'^i^LL)ij i^^'^ similarly for the RR 
sector), {mlLR)^^ = Addiag{md,ms,mb) and (m^ j-^)^^. = m1{S'[j^).. = m|((5^^)*., where 
i,j = 1,2,3 and i ^ j. (Analogous parametrization is used also in the up-sector.) In the 
case of 6 — > 57, the relevant mixings are those between the 2nd and 3rd generation squarks. 
For simplicity, in this section we will use the notation = {^'[1)2^ ete. In all the cases 
presented in this section, we keep only one 5^'" non-zero at a time and set all other (5^'^'s 
to zero. 
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Figure 4: The Wilson coefficient Cj{iiw) vs. tan/3 for A/susy = 1 TeV, rrig — \/2mq ~ AfsusY, 
fi = 0.5 TeV, Ad ^0, Au = -ITcV, mH+ = 0.5 TeV and Sf^ = 0.2. All the other S'^-'^'s are set 
to zero. In the upper left (right) window we show the chargino (gluino) contribution to C7(/iw), 
and its total value from new physics in the main (lower) window. The dashed line shows the case 
of LO matching, the dot-dashed one includes the effect of NLO QCD, while the solid one includes 
also the effect of beyond~LO corrections from squark-gluino loops. 



First, we present the effect of including beyond-LO corrections derived in the previous 
sections. In the three panels of fig. ^ we plot the Wilson coefficient C7(//w) vs. tan/3 for 
Sf^ = 0.2, fj, = 500 GeV and for other relevant parameters as specified in the figure caption. 
In particular, is chosen such that an approximate cancelation between the charged Higgs 
and the chargino contributions takes place in the case of fi > for = 0, similarly to 
the case of the Constrained MSSM (CMSSM). The dashed (solid) line corresponds to the 
LO (beyond~LO) value. (Here the LO matching is defined so that the SUSY contributions 
8 are evaluated at //W) and the NLO corrections S^-Cj^l and ^Y'^Cj^g are neglected.) 
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Figure 5: The Wilson coefficient Cjipw) vs. tan/3 for (5^^ = 0.05 and for tlie otfier parameters 
as in fig. ^. 



The dash-dot hne indicates the effect due to the NLO QCD correction only. Relative to 
the LO approximation, the chargino contribution, which in this case dominates Ci{fiY/), is 
significantly reduced towards zero. As a result, the overall value of C7{fi-w) is also reduced, 
although, like in the LO case, it does grow with tan /3. (One finds a fairly similar effect for 
Csl/^w) as weh.) 

Next, in fig. ^ we plot (^^(/iw) ^s- tan/3 for dj^i^ = 0.05, /i = 500 GeV and other pa- 
rameters as in the previous figure. In this case, at LO the chargino contribution is zero 
while beyond-LO effects generate it on the negative side, with the magnitude growing with 
tan (3. In this sense, like before, the beyond-LO chargino contribution is again significantly 
reduced relative to the LO value. However, at some point it out-balances the gluino con- 
tribution which is also reduced but remains on the positive side and is roughly independent 
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Figure 6: The Wilson coefficients C7{^y^) (upper left) and Csi^J-w) (upper right) vs. tan/3 for 

S'rl = 0.' 
as in fig. 



.05; and Cy(/iw) (lower left) and Cg(/xw) (lower right) for Sf^j^ — 0.5. All other parameters 



of tan/3. The total C'-j{^^n) in this case decreases and at some point even changes sign. On 
the other hand, the total Cg(//-w) is reduced by a factor of about two but remains positive 
and roughly independent of tan/?. The main effect in Cg(//w) comes from reducing the 
gluino part, while the chargino one is only slightly reduced down from zero. 

Substantial beyond-LO effects also appear in the LR and RR sectors. In order to 
illustrate this, in the four panels of fig. ^ we plot Cr^sifJ'w) and C7g(/xw) vs. tan/3 for 
^LR ~ 0-05 and Sf^j^ = 0.5, respectively. As a rule, Csi^w) and Cg(//w); like in the previous 
cases, are reduced mostly because of a strong suppression of the gluino contribution. In 
all the cases of Cj^sifJ-w) and presented in fig. ^, one finds that, in addition, the 

chargino contribution shifts from a positive, or zero, value down to an increasingly negative 
one, although much less so for Csi^J-w) than for Cj{ij,t^). 

In the case of /i < 0, the two competing effects from the RG running and from the 
gluino-squark corrections lead to the NLO-corrected Wilson coefficients either decreasing 
or increasing with tan /?. This can be seen in the four panels of fig. where all the other 
parameters are kept as before. In general, one can see that beyond-LO corrections to the 
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Figure 7: The same as in fig. I but for fi < 0. 



Wilson coefficients are substantiaL However, in contrast to the case of > 0, the Wilson 
coefficients can now either be larger or smaller relative to the LO approximation. 

For comparison of the general GFM scenario with the limit of MFV, in fig. ^ we plot 
CrilJ-w) ^^s- tan/3 for the case of MFV and for n = ±500 GeV and the other parameters 
as before. The coefficient Cs{fiw) shows a similar behavior. It is clear that the effect of 
including beyond-LO corrections is not as striking as in GFM. 

Large corrections to the Wilson coefficients in the case of GFM (but not MFV) often 
leads to substantial changes in the predictions for BR{B — > ^57). In fact, the overall 
tendency is to reduce the magnitude of super symmetric contributions relative to the LO, 
especially for large tan/3 and /i > 0. This leads to the "focusing effect" of concentrating 
on the SM value which we have identified in the previous paper [26|. We illustrated the 
effect in fig. ^ where we plot BR{B Xs^) vs. rriq for some typical choices of parameters. 
The focusing is indeed rather strong. In contrast, in the case of MFV the effect is much 
less pronounced p^ . 

As explained in ref. [26| and in sec. |l|, focusing originates from two sources. Firstly, 
renormalization group evolution of Cy^s from /UsusY to /iw reduces the overall amplitude 
of these coefficients. In particular, the reduction of the strong coupling constant at /xgusY 
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Figure 8: The coefficient C7(/^w) vs. tan/3 in MFV for < (left) and ^ > (right). All other 
parameters as in fig. ^. Relative to GFM, the difference between LO and beyond-LO effects is 
rather mild. 



has substantial effect on the gluino contribution. The NLO QCD matching at /^SUSY 
brings further suppression. Secondly, the gluino loop contribution to the mass matrix of 
down-type quark and the leading NLO SUSY-QCD corrections to the Wilson coefficients 



5f^ci'l^^^ have a common origin. This considerably reduces the LO gluino 



-"7,8 

contribution to Cjg, depending on the sign of /x. At ^ > this correlation works as 
alignment, which is enhanced by tan /? and explains the bulk of the focusing effect. For 
/i < 0, SUSY-QCD corrections cause anti-alignment instead, which competes with the 
overall suppression by renormalization group evolution and results in small focusing (or 
even de— focusing) . Some NLO suppression of SUSY contribution already exits in MFV. 
However, flavor mixing is essential for the focusing effect with GFM as described above. 

One important consequence of the focusing effect is a significant relaxation of bounds 
on the allowed amount of mixing in the squark sector relative to the LO approximation. 
We illustrate this in fig. f where we delineate the regions of the plane spanned by df^^ and 
niq where the predicted values of BR{B Xs"f) are consistent with experiment, eq. (p..lj), 
at Ifj (long-dashed) and 2a (solid). One can see a significant enlargement of the allowed 
range of relative to the LO approximation. Similar strong relaxation due to NLO-level 
corrections is also present for 5'^^^ and (5^j:j, and somewhat less so in the case of d'^j^. 

One should also note that, with even a small departure from the MFV case one can 
significantly weaken the lower bounds on rriq ~ MgusY derived in the limit of MFV. 
The effect is already present in the LO approximation but it generally becomes further 
strengthened by including the beyond-LO corrections considered in this paper. It results 
mostly from a competition between the chargino and gluino contributions, the latter of 
which is in MFV absent. 

The fact that lower limits on MgusY obtained in the case of MFV are highly unstable 
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Figure 9: The branching ratio BR{B Xs7)/10 ^ vs. rriq 



-^SUSY for in the LO and 



beyond-LO approximation for Sfj^ = ±0.2 (upper left), (5^^ = ±0.05 (upper right), (5^^ = 0.05 
(lower left) and (5^^ — 0.5 (lower right). All other parameters as in fig. |4[ 



with respect to even small perturbations in the defining assumptions (2.34) of MFV, ap- 
pears to be rather generic. To show this, in the four panels of fig. |ll| we plot contours of 
BR{B Xg^) in the planes of 5^ and rriq for tan /3 = 40, /i > and all the other param- 
eters as specified in the figure caption. The bands delineated by dashed (solid) lines mark 
the la (2(t) regions consistent with experiment ( |1.1| ). As before, only one 6'^''^ is kept dif- 
ferent from zero in each case. Note that in the regions in between the two allowed bands in 
the case of df^ and df^^ the branching ratio reaches a shallow minimum whose position de- 
pends on the relative cancellation of the chargino and gluino contributions, and which also 
depends on tan /? and other parameters. A similar shallow minimum appears in the case 
of small df^^ or and smaller rriq. (It is worth noting that the values of BR{B Xg'^) 
in the excluded region between the two bands is are such that a fair decrease of the ex- 
perimental value would considerably reduce it.) Note that the cases of (5^^ and the 
plots are almost symmetric. This is because of the new gluino contributions that appear in 
Cy g and which contribute to the branching ratio quadratically without interfering with the 
MFV contributions in Cy^s- In these cases, new contributions always come constructively. 
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Figure 10: Contours of the branching ratio BR{B Xs7)/10^"' in the plane of rriq — mgq — 
■^SUSY and S'^j^ in LO approximation (left) and with dominant beyond-LO corrections included 
(right). Regions of the plane of rriq and (5^^ where the branching ratio BR{B Xs"i) is consistent 
with experiment, eq. ( |l . l| ) , at Icr (2(t) are delineated with a long-dashed (solid) lines. We take 
tan/3 = 50, mg = \/2mq, Ad = Q, Au = —niq, /i = 0.5 m^, and m//+ = 0.5 m^. The magenta 
short-dashed line marks the case of MFV. All the other J'^'s are set to zero. 



For fi < (not shown in the figure) almost all of the region is excluded. 

In the plots, BR{B Xg'^) appears to be more sensitive to and (5^^ than to 
d'^i^ and in our phenomenological approach. However, note that, once we assume a 
correlation between the A-terms and the Yukawa couplings, like in the MFV, the natural 
scale of (<^^l) reduces from one to mfe/MsusY < 0(10-2 )_ 

Even though the focusing effect considerably reduces the gluino contribution in the 
case of GFM, the branching ratio still shows strong dependence on Sf^j^ and df^^ because 
of the interference between the gluino contributions and the MFV contributions to Cy^g- 
In addition, a new LO chargino contribution, other than the ones from the CKM mixing, 
is present in the case because of the SU (2) l relation. 

It is often claimed that the case of /i < is inconsistent with BR{B ^sl)expt (Li)- 
This is so because in the case of MFV the chargino contribution adds to the SM/2IIDM 
contribution constructively for ^A^ > and in unified models, like the CMSSM, with the 
boundary conditions Aq = (or similar), the RG running generates Au < at mw- As a 
result, theoretical predictions are inconsistent with experiment, unless MsusY ^ 2TeV, or 
so. 

However, the claim is only true in the case of the strict MFV but not necessarily for 
even small departures from it, as can be seen in fig. 12. Indeed, for ~ —0.1 and/or 



(5^^ ~ —0.02 one can even evade the bounds basically altogether. The effect is due to a 
partial cancellation of the chargino and gluino contributions which in the case of 5^^, can 
be efficient enough already in the LO approximation but in the case of df^^ would only 
allow rriq down to some 500 GeV but not lower. In contrast, beyond LO one can have 
as small as some 200 GeV, as can be seen in fig. 12. On the other hand, for no ranges of 
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Figure 11: Contours of the branching ratio BR{B Xs7)/10~'* in the plane of rriq = mgq = 
-MsuSY and <5f ^ (upper left), Sfj^ (upper right), Sf^^ (lower left) and Sj^j^ (lower right) for tan/3 = 40 
and rrig = rriq, Ad = 0. All the other parameters are kept as in fi g, p^ . The long-dash (solid) curves 
delineate the regions which arc consistent with experiment, eq. (1.1), at Icr (2(t). 



—0.15 < dj^j^ < 0.15 or —0.8 < < 0.8 and nig < 1.5 TeV could we find any region 
consistent with experiment at either LO or beyond-LO level. 

In tables we present allowed ranges of the off-diagonal entries for a number of 
cases. As can be seen in figs. 1^12, typically there are two bands. In the cases of S'li^ 
and (5^^, one is around zero and one increasingly deviating from zero as -A^fgugy increases. 
On the other hand, the allowed bands of (5^^^^ and Sf^^ tend to be roughly symmetric, as 
mentioned above. As one can see from the tables, the allowed bands do not necessarily 
increase with tan/3. This is because the contributions from the chargino and the gluino 
show a different dependence on tan/3, and their approximate cancellation, can happen at 
different values of tan /?, as for example fig. |5| illustrates. 



The significant enlargement of the allowed parameter space in the framework of GFM, 
relative to MFV, leads to a strong relaxation of constraints on the allowed ranges of the 
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Figure 12: The same as in the upper two windows of fig. p]but for /i < 0. In contrast, no ranges 
of (5^^ and (5^^ (lower two windows of fig. [ll| ) are consistent with experiment. 
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Table 1: Allowed range of (5^^ beyond LO, which resides within the 2a experimental error. Pa- 
rameters other than specified in the table are fixed to Mi = M2 = rriq, mt — 0.5mg, Au — —rriq 
and Ad = 0. 



"bare" CKM matrix element K^^^ relative to K^b — 0.04 which is determined by experi- 
ment. In the two windows of fig. 13 we present contours of the ratio / Kch in the plane 
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Table 2: Allowed range of (5^^ with the same parameter set as in Table |l|. 
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Table 3: Allowed range of (5^^ with the same paramter set as in Table H. 



fj, = 0.5 rrig, and mff+ = 0.5 rUg and show two cases tan/3 = 40 (left) and tan/3 = 60 
(right). One can see that the ratio /Kch can exceed a factor of ten, or more, for large 
enough, but still allowed, ranges of (J^/j- It is also interesting that the ratio can be much 
smaller than one and that it can even "cross" zero^^ in the allowed region between the two 
contours of 1. In this case the measured value of is of purely radiative origin and is 
generated by GEM in supersymmetry On the other hand, for large < 0.8 the 

ratio does not exceed a factor of three. 



By definition, if^j is positive in the standard notation [}43|. Therefore the sign change can always be 
rotated away by redefinition of (s)quark fields and does not appear in the figure. 
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Figure 13: Contours of the ratio Kf^ jKch in the plane spanned by m-g = mgq and (5£^. We fix 
rrvg = ^/2mq, Ad = 0, A.^ = —mq, /i — 0.5 mg, and mff+ — 0.5 and show two cases tan/3 = 40 
(left) and tan/3 = 60 (right). 



6. Conclusions 

The squark-gluino loop correction that appears beyond the LO in the inclusive process b 
s'y shows a large effect in the MSSM with general flavor mixing. Its main feature is that of 
reducing the magnitude of supersymmetric (mostly gluino) contribution to BR{B — > ^^7) 
relative to the LO approximation. This focusing effect leads to a considerable relaxation 
of experimental constraints on flavor mixing terms and on the lower bounds on MsusY 
obtained in the limit of MFV, to the extent of even allowing small MsusY at /i < 0. 

Given such a large effect appearing beyond LO, one may question the validity of 
perturbation theory in the analysis presented here. However, in this case the mechanism 
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that generates it, namely the large gluino-squark corrections, is absent at LO. We believe 
that the mechanism is a dominant one at NLO, although we have examined it within the 
limitations of our assumptions. It would therefore be highly desirable to have available a 
complete NLO calculation of BR{B — > Xs^) in supersymmetry with GFM which would be 
applicable in more general circumstances than those considered here. It is also clear that, 
with improving experimental precision, an NLO-level analysis appears to be essential in 
testing SUSY models of flavor at large tan/3, like, for example, those based on a SO(IO) 
SUSY GUT Pi. 
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A. Effective and Vertices 

In this Appendix we summarize the effective vertices of PF^-boson (unphysical scalar G^) 
and charged Higgs boson which arise from integrating out gluino-squark loops at the 
matching scale fisuSY- All the coupling constants and masses appearing in the following 
are evaluated in the effective SM at some arbitrary scale n although in actual expressions 
for Wilson coefficients they will be evaluated at the scale usuSY- We work in the MS- 
scheme except for the case explicitly mentioned. The renormalization group equations for 
evolving the couplings in the effective SM to ^isusY from their values at are given 
in [^]. We work in the physical super-CKM basis where the loop-corrected mass matrices 
"f^d^u (|2.3|)-(^7^ of both the down- and up~type quarks are diagonal. 



A.l VF='= boson 

The effective udW vertex, evaluated at some scale /x, is given by 

3 

^ = E U^{^Ji)l^W,PLKtfd, + h.C. (A.l) 

jj'=l 

(A.2) 

Note the effective CKM matrix KfJ which appears in ( A.l ), instead of Kij = 1,2,3) 



K'rf = K,, + f ^ ) AK,„ (A.3) 



where 

3 6 

AK,, = 4C72(3) J2 E iKL)H{^uL)ikKH{nL)ji{^dL)j,C24{ml,ml,m'~^;fi) 
k,i=ii,j=i 

1 ^ 

+ 2 E [(^KL)HKkj + K,UAZ,L)kj\ ■ (A.4) 

k=l 

The first term above is due to the VT-boson coupling to the intermediate squark fields, 
while in the remaining two the VF-boson couples to the external quark lines, and 

6 

{AZ^Lhjifi) = 2C2{3)Y,iKL)H{'^uL)IJB^{m|,ml,0■,^^), (A.5) 

1=1 

6 

{AZ.Lhif^) = 2C2(3)^(rSJ,,(r,i),^.i?i(m|,m|^,0;M). (A.6) 

1=1 

It is the elements Kfj^ that are equal to the experimentally determined elements of the 
CKM matrix. The elements of Kij are then determined by ( [A.3| ) but the difference with 
Kfj^ is rather tiny. 

Note that the effective vertex defined by ( |A.1| ) gives finite contributions since In /i and 
infinite terms which appear in the Passarino-Veltman functions, which are collected in 
appendix 0, disappear in actual calculations. 
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A. 2 Unphysical scalar 

The effective vertices of the unphysical scalar with d and u, evaluated at the scale fiw 
in the 't Hooft-Feynman gauge, are given by 



£ = ^ (cr" ^^'l + C^^' ^Pr) G+dj + h.c, 



(A.7) 



where the couplings (f^"^^^ receive contributions from the tree-level q'^"-^^^^^ from 



squark-gluino loops (j^^^'^^^ 



(Ai 



The former are given by 



^Gu,d,{0) ^ g2Kijmdj{^Xy^) 



\/2mw ' \/2mw ' 

where rnujinw) and fridjinw) denote the running quark masses at /i^y- For simplicity we 
have neglected here a small rotation of the quark mass basis due to the renormalization 
group evolution between and ;^SUSY- 
The corrections are given by 



^Guidj{l) 



^Guidj{l) 



3 

r/7 



/,J=i 

3 



3 

r/7 



k,l,m=l 
6 



(A.IO) 



i,j=i 



2^2(3) {r:L)n{rdR)jAc'''-'~^') mgCo(m|^,m| m|) 



k,Lm=l 



(All) 



where Cq is a Passarino-Veltman function and the expressions for {5m(i^u)ik are given in 
and below eq. ( |2.26| ). 
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The factor 77 = o;cj(//susY)/as(/^w) appears as the result of running the couphngs 
^LR*^^^^^' which are initially evaluated at nsvsY, down to /iW; with the SM QCD RGE's 
including six quark flavors. 



The functions AZuL and AZ^i are defined in (A. 5) and ([A.6D, and 



92 



+ 



92 



Ijm 



k,l,m=l 
3 



-7=^ X] (^dL)ikKlk 



u,LR)i^ - <'<'ui 



niuM cot f5 Sim 



(r:ii)j^.(A.12) 



Like in the case of the effective coupling of the M^-boson above, In and infinite terms 
in Bq effectively do not contribute because of the unitarity of the diagonalization matrices 



A. 3 Charged Higgs boson 



The effective vertices with d and u, evaluated at the iJ,-yv, are given by 



Huidj 



(A.13) 



where, similarly to the case of the unphysical scalar before, the couplings c ' receive 
contributions from the tree-level (j^^'^i^'^^ g^j^d from squark-gluino loops (j^^'^i^^^ 



^L,R + 4^ '~'L,R 



(A.14) 



The former are given by 



C 



V2mw 



C^u.d,io) ^9_2K^^^^^^ 



'R 



V^mw 



(A.15) 
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where rriui and m^i have been defined below eq. (|A.9| ). The corrections are given by 



Huidj{l) _ 3 



l,J=i 

02 cot ^ A 



(A.16) 



Huidj{l) _ 3 



R 



l,J=i 
3 



(72 tan /3 



+ 



2 2\ 

m- ,my.) 



92 "'an /J ST^ , . *n 



fc,Z,m=l 



(A.17) 



As in the previous subsection, the factor r] appears as the result of running the couplings 
C^^^^^^\ which are initially evaluated at /^susY, down to /iW; with SM QCD RGB's 
assuming six flavors. 



(jHdiuj 



92 



92 



k,l,m=l 
3 



Jm 



+ 



-j^ — {TdL)ikKik ("^^,L/^);^cot/3-/im„;5zm (r^^) j^. (A.18) 

k,l,m=l 



B. Effective x^i X° g Vertices 

We summarize here the effective vertices of the chargino neutralino ^iicl gluino g 
fields, evaluated at /UsusY- Notice that the corrections to the tree-level vertices come as 
a result of using gluino-squark loop-corrected values for the quark and squark quantities 
appearing in the expressions below. 

B.l Chargino 

The effective quark-squark-chargino vertices are given by 

3 2 6 

^ = Y.Y.Y.^^^^^'iL^ialPL + {CdR),alPR\Xa^I + h.C., (B.l) 

i=l a=l 1=1 
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where 



3 C^Wn 



\/2mH/ COS /3 



92 (m*-"''*) 



(B.2) 
(B.3) 



and V and C/ are the 2x2 two matrices with which one diagonahzes the chargino mass 
matrix. 

B.2 Neutralino x° 

The corrected quark-squark-neutrahno vertices are given by 



3 4 



(B.4) 



i=l r=l 1=1 



where 



{NdR)irI 



-V2g2 



-tan.^A7,(r^,), + A_^^Pr3(r*,),, 

iiV^s + ^tan^H-iVrii) {nL)ii 



(B.5) 



+ 



/ (0)*N 



ki 



2mw cos P 

and A'^ is the 4x4 matrix with which one diagonahzes the neutrahno mass matrix |34| 
B.3 Gluino g 

The effective quark-squark-gluino vertices are given by 

3 6 



(B.6) 



= EE*^^ i(GdL)uPL + {GdR)uPR]dig 

i=l 1=1 
3 6 

+ [{Gul)uPl + {Gur)uPr] uig + h.c, 



i=l 1=1 



where 



dLJliJ 



(B.7) 



(B. 



and analogously for the up-type sector. 

The above effective couplings for the chargino, neutralino and gluino reduce to their 
tree-level values by making the following replacements: {TduLR)ii — *■ (j'du.LR) ^^'^ 



m 



(0)^ 
d,u , 
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C. Passarino-Veltman functions 



We present here explicit expressions for the Passarino-Veltman functions [45| which appear 
in our calculation. We take all the external momenta equal to zero. 



Co{x,y,z) 
C24{x,y,z;fi) 



1 xlnx-ylny 

— him/i — 1 

e X — y 

11 , 13x-y 

mu — 7 

2 e ^ 4 x-y 

^2 



+ 



X- 1 (2x - y)y 

in X — — — in y, 



{x - yY 

x\nx 



2 {x-yY 



+ 



{y - x){z - x) 
y\ay 



+ 



zlnz 



{z-y){x-y) {x-z){y-z) 

Inx 



1 r 1 , 3 

4|- + 21n^+- 
Iny 



(x - y){x - z) 
z^ In z 



{y-x){y-z) {z-x){z-y) 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



D. Mass functions for the Wilson coefficients 

In this section we collect several auxiliary functions which appear in sec. g The mass 
functions for the SM and 2HDM contributions to Cy^g are given in ||8|, 



(1) , ^ _x(7 - 5x - 8x^) ^ x^(3x 



24(x - 1)3 



+ 



4(x - 1)4 
3x2 



■ Inx, 



(1) . . _x(2 + 5x - x'^) 
8 8(x - 1)3 4(x - 1)4 



Inx, 



Fi'^\x) =I^^JL^ _^ x(3x - 2) 



Inx, 



12(x-l)2 6(x-l)3 

(2) _ x{3-x) X 

8 ^'^^-4(x-l)2 2(x- 1)3 



(D.l) 
(D.2) 
(D.3) 
(D.4) 



The mass functions and related constants for the x ) X ^iid g contributions to C7, 
are given in [41|, 
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iff! 



-"3 



HP 



[7] 



-^6 



InxH ^^^^t:; ^2 — , (U.5) 



6(1 36(1 -x; 



X -x^ + 5x + 2 

InxH — — (D.6) 



2(1 -x)4 12(1 -x; 

-3x2 _^ 2x , -5x2 

lnx + — (D.7) 



^ 3(1 -x)3 6(1 -x)2 ' 

X — x2 + 3x /T^ X 

")=(r3^^""+2(r3^' (^-'^ 

x)=-^Ff](x), (D.9 

x)=Hf\x), (D.IO) 

x)=-U//fl(x) + ^V (D.ll) 



x) =/7^^'(x) + -, (D.12) 



2' 

-iFf](x), (D.13) 

9x2 -X 19x2 + 40x-ll 

■InxH z -1, , (D.14) 



16(1 -x)'' 96(1 -x)3 

:-i(i7P(x) + i), (D.15) 

Inx + i^^, (D.16) 



' 8(1 -x)3 8(1 -x)2' 

Al^l =-, At^l = i. (D.17 
6' 



E. NLO QCD corrections to the Wilson coefficients 



For completeness, in this section we collect the NLO QCD corrections to the Wilson coef- 
ficients given in js], 41]. 

The NLO QCD corrections to the Wilson coefficients in the SM are summarized as B, 



15 + 6 In 
0, 

V "^w J i i 



+ 



m?(Mw) 



In-^. 



for 
for 
for 
for 



i = 2, 3, 5, 6, 



4, 

7,8, 



(E.l) 
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where 



A7(x) 



A8(x) 



^(2;)=^ 7 ^ 7 71 ^Inx Inx, (E.2) 

^ ^ 12(x-l)3 6(x-l)4 3 ' ^ ^ 

-lex'^ -122x3 + 80x2 -8x^. / I 

= 90^31? V x 

6x^ + 46x3 - 28x2 2 

H 7 TE In X 

- 1)^ 3) 
-102x5 - 588x'^ - 2262x3 + 3244x2 - 1364x + 208 , 

H ; TE Inx 

81(x-l)5 

^ 1646x'^ + 12205x3 - 10740x2 + 2509x - 436 



486(x - 1)4 
208 - mix + 1086x2 + 383x3 _^ ^2x'^ 



81(x - 1) 



4 



(E.4) 

2x2(14 -23x- 3x2) , 
+ 3(^31)5 1-^' 

^ , , -4x4 + 40x3 + 41x2 + x^. / 

= 6(^^T? - X 

-17x3 - 31x2 , 2 

H ; — in X 

2(x - 1)5 

-210x5 + 1086x4 + 4893x3 + 2857x2 - 1994x + 280 , 

H ; Inx 

216(x - 1)5 

737x4 _ 14102x3 - 28209x2 + 610x - 508 
+ 



1296(x - 1)4 
140 - 902x - 1509x2 - 398x3 _^ 773.4 

~ (E.6) 

x2(31 + 17x) , 

H ^7 r^lnx. 

2(x - 1)5 



Similarly, the NLO QCD corrections to the charged Higgs boson contributions to the 
Wilson coefficients are given in 



5fcf)(Mw) 



0, 

E" {y) , 

Gf (y) + Af (y)lni^ 
Gf (y) + Af (y)ln 



(y) 

\E^ (y) 



for 
for 
for 
for 



1,2,3,5,6, 

4, 

7, 



(E.7) 
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where, y = m'^{fiw)/mjj, 



A?(y) 



Gf{y) 



yil6 - 29y + 7y^) , y(3y - 2) 



36(y - 1)3 



+ 



6(y - 1) 



■Iny, 



2 

+ -y 

9^ 



4(-3 + 7y- 2y^ 
3{y - 1)3 

2(-3-y + lV 
3(y-l)4 

y(18-37y + 8y2) 



+ 



+ 



-50 + 251y - 174?/ 



1\ 8-14y-3y2, 2 

2y3) 7_i3y + 2y2 

(y-l)3 
y(-14 + 23y + 3y2) 



1-i 

y 

2 im.,3 



192y^ + 21y^ 



9(y-l)5 
797 - 5436y + 7569y2 - 1202?/3 



(y - ir 

Iny 



108(y - 1)4 



2 

+ 7Ty 



3^ 



1 

+ 72/ 
b 



21-47y + 8y2 2(-8 + 14y + 3y2) 

(y-l)3 + ^^^^ 

-31 - 18?/ + 135y2 



(y - 1)^ 

14y3 y(14 - 23?/ - 3?/2) 



6(2/ -1)4 
-36 + 25y - 17y2 



+ 



2(y-l)3 
-3 - 187y + 12y2 



4(y - 1)4 
y(30 - 17y + 13y2^ 



+ 



i-i 

y 

Iny + 



Iny 



{y - 1)5 
y(19 + 17y). 2 
^^^^"^ 
3(143 -44y-F29y^)" 
8(y-l)3 . 



y(31 + 17y) 2 
In^y 



+ 



+ 



(y-l)4 yj (y-l)5 

226 + 817y + 1353y2 318y3 42y4 



36(y - 1)5 
1130 - 18153y + 7650y2 - 4451y3 



■Iny 



216(y-l)4 



;y 



1 

+ 

6 



81 - 16y + 7y^ 

2(y-l)3 
-38 - 261y + 18y^ 

6(y - ir 



19 + 17y ^ 
^Iny 



(2/-I) 



7y3 y(31 + 17y) 



(Ei 



In^y (E.9) 



(E.IO) 



(E.ll) 



(E.12) 



The NLO matching conditions for 6 ^ 57 in the MSSM have been calculated in the 
limit of nig — > 00 by Bobeth, et al, in [41|. We apply their calculation on the 2-loop gluon 
corrections to Cy^s at /^SUSY; while distinguishing between the two scales //w and nsvsY- 
We neglect the subdominant contribution from C4 and do not include the corrections from 
quartic squark vertex in light of the different assumptions about the mass hierarchy. 
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The NLO QCD correction to the chargino contribution to is given in |41|, 



C'?8(/^SUSy) 



1 ^ ^ 



m 



w 



m 



Xa 



{CdR)2al{CdR) 



3al \ ^1 



H 



[7,8]'/ u! 



Xa ^ 



[7,8]", 



x^'_ ) In 

■ Xa 



■ -^{CdR)2al{CdL)laI 

mi) 

x|Ff'«l'(x"^_)+/?p]"(x"^)h. 



( ^SUSY A 



(E.13) 



where x^'_ = muj"^ / m^-'^ and 



Xa 



hP"{x) 



24x3 + 52x2 - 32x, 



Lis 1 - - + 



-189x3 - 783x2 _^ 4253. _^ 43 



+ 



9(1 - x)4 
-1030x3 - 1899x2 + 1332x + 85 
243(1 -x)4 ' 



81(1 -x) 



Inx 



6x3 _ 523,2 _^ 322; 28x3 - 129x2 - 12x + 41 

: ; mX H 

9(1 -x)5 
112x2 - 48x, 



9(1 - x) 



Li2 1 - - + 



27(1 -x)4 
1\ 12x3 - 176x2 + 64x + 16 



+ 



-170x2 + 66x + 20 



9(1 - x)3 
12x3 - 88x2 _^ 4Q^ 
9(1 - x)4 

-9x3 _ 4g^2 _ 4g^ 



Inx + 



9(1 - x)4 

14x2 - 54x + 32 
9(1 -x)3 ' 



Inx 



+ 



f l\ 81x3 + 594x2 + 1270x + 71 

12(1 -x)4 V ~ ^ y 108(1 -xf 

923x3 + 3042x2 + 6921x 1210 



Inx 



648(1 - xY 



5x2 _^ 7^3 _ 303,2 _^ ^4^3. _^ 26 
Inx -1 



3(1 -x) 



18(1 - x)4 



(E.14) 
(E.15) 
(E.16) 
(E.17) 
(E.18) 
(E.19) 



[8]/, , -16x2 -12x^. 1\ , 52x2 + 109x + 7 



+ 



95x^_+180x+_61 
12(1 -x)3 ' 



(E.20) 



[8],, 10x2 + 26x ^2 + 30^ 

^2 (")= 3(l-x)4 ^""+ 3(1 -x)3 • ^^-''^ 



^^Strictly speaking, higher-order SQCD correction is also included in the formulae below through gluino- 
squark corrections to the squark vertices and masses. 
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The NLO QCD correction to the neutrahno contribution to C7.8 is given in |41] 



6 4 



+ ^{NdR)2rl{NdL) 



m 



where x^n = "t-j ^/m^o^ and 



(E.22) 



16x2 + 282;^ / 1\ -108x2 - 358x - 38 , 

: — -, 71— L12 1 H ; TP Inx 

9(1 -x)4 81(1 -x)5 



+ 



23x^ - 765x2 _ g93^ _ 

243(1 -x)4 ' 



4x2 - 28x , -8x3 _^ 42^2 _ §4^ _ 22 
Inx H 



9(1 -x)5 
16x2 + 48x 
■ 9(1 -x)3 



27(1 -x)4 
^ , 1\ -8x2-68x-8, 
^^^'^--'+ 9(l-x)^ 



+ 



-26x2 - 54x - 4 
9(1 -x)3 ' 



^4 ^"^^ 9(1 -x)4 
i/0(x) =i/g'(x), 
<'(x) =<'(x), 



8x2 _ 443, ;^Q^2 _ _ -^^ 

Inx + 



9(1 -x) 



(E.23) 
(E.24) 
(E.25) 

(E.26) 

(E.27) 
(E.28) 



The NLO QCD corrections to the gluino contribution to C7 § is given by 



-^KsVsusy) 



w 



H {GdR)2l{GdL)3I 

rub 



(E.29) 
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where = ^/m| and 
g di I 9 

„\7],, ^ 17x2 + 86x- 15^. / 1\ 6x3 + 45a;2 + 66x-5, 2 

HI ^ (x) = -; Ll2 1 H -, m X 

^ ^ ' 18(1 -.t)4 ^\ xj 12(1 



+ 



-36x'^ - 315x'^ + 1161x2 + 751x + 23 
162(1 - x) 

'799:r^ + 1719.C- + 10131.r - 1847 



+ ttt^t; — In X (E.30) 



972(1 - x) 



„n\n, , 18x3 + 107x2 + 43x , -Sx^ + 384x2 + 699x + 20 

^(13^)5 1"^ + 108(1 - xY ' (^-^'^ 

„r7i/, , 19x2 + 60x - 15^ / 1 \ 6x3 + 35^2 _^ 43^ _ 5 

-"fi (a^) = TTT^ ^ L12 1 H TT- ri m X 

^ ^ 9(1 -x)3 ^ V x/ 6(1 -x)4 



-27x3 + 106x2 + 52x + 1 14x2 + 333x - 83 

+ 31 Inx + ,0/. _^^3 ^' 



9(1 -x)4 18(1 -x) 

mi, , 18x3 + 80x2 + 28x , 55x2 + 69x + 2 

— InxH ; 

9(1 - x)4 9(1 - X 

[8]/_ 45x3 - 1208x2 + 901x- 570 _ 1 



(E.32) 



<'"(x) = Zrj, - 1,^+ — J— 7 , (E.33) 



'H^ 

-237x3 - 846x2 + 282x - 95 , o 

+ '"^ 

2520x^ - 10755x3 - 10638x2 - 6427x - 44 

^ 864(1 - x)5 

5359x3 _ 241425x2 + 143253x - 59251 

^ 5184(1 -x)4 ' 

..[81/// X -747x3 - 640x2 + 43x , -779x3 _ 7203x2 - 93x + 11 

Hl^ (x) = TT mxH -, —. , (111. 35) 

^ ^ ' 48(1 -x)5 288(1 -x)4 ' ^ ' 



^[8]/, , -359x2 + 339x- 204^. / ^ 

= 24(13^^ V X 

-78x3 _ 333^2 _^ 2Q5^ - 34 , 2 

H ; 71 In X 

8 1 - X 4 

„ o (E.36 

-207x3 - 1777x2 + 23x - 151 , 

H ; n inx 

48(1 - x)4 

^ -1667x2 +990x- 379 



24(1 - x)3 
- 133x2 + 
6(1 - x)4 "^"^ ' 12(1 -x)3 



rg],, -126x3 - 133x2 + 7x. , -553x2 + 84x - 35 
^6 (^) = FTi ^4 lnx + — -3 . (E.37) 



Similar expressions can be obtained by interchanging the indices L and R 

in the above formulae. 
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